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SUMMARY 

The paper presents the results of a numerical investigation of the effects 
of buoyancy on the direct condensation of saturated or nearly saturated vapor 
on a horizontal liquid surface in a cylindrical tank. The liquid motion 
beneath the liquid-vapor interface is induced by an axi symmetric laminar Jet of 
subcooled liquid. Analysis and numerical results show that the dominant param- 
eter which determines the influence of buoyancy on the condensation rate is the 
Richardson number. However, the effect of buoyancy on the condensation rate 
cannot be quantified in terms of Richardson number alone. The critical value 
of the Richardson number below which the condensation rate is not significantly 
reduced depends on the Reynolds number as well as the Prandtl number. 
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NOMENCLATURE 


area, jet, outflow or interface 
tank to jet diameter ratio, D/d 
specific heat-at-cons tan t-pres sure 
tank diameter 
jet diameter 

Grashof number, g|3(T s - T^)D^/v 2 
acceleration due to gravity 

local condensation heat transfer coefficient, m c h^g/(T s 
average condensation heat transfer coefficient, n^h^ /(T 
latent heat of condensation 
Jakob number, Cp(T s - Tj)/h^.g 


s 



V 


‘Member, AIAA. 


k thermal conductivity 

m c ,m c local and average condensation mass flux 

p pressure 

Pg equilibrium hydrostatic pressure 

2 

p* dimensionless pressure (p - p )/pu, 

3 W 

Pr Prandtl number, yC /k 

r 

r radial coordinate 

r* dimensionless radial coordinate, r/D 
Re^ jet Reynolds number, pu^d/y 

2 2 2 

Ri Richardson number, g|3(T s - T\)D/Uj or Gr/B Re^ 

St c local condensation Stanton number, ^/pu^Cp 
St average condensation Stanton number, h /pu.C 

C C J p 

T temperature 

T* dimensionless temperature (T - T.)/(T - T.) 

J S J 

u axial velocity 

u c condensation Induced velocity 

u*,u* dimensionless velocity, u/u., u /u. 

C J C J 

v radial velocity 

v* dimensionless radial velocity, v/u. 

J 

v* dimensionl ess- radi al velocity at the Interface, v /u. 
* ^ J 

x axial coordinate measured from the tank bottom 

x* dimensionless axial coordinate, x/d 

3 coefficient of thermal expansion 

y dynamic viscosity 

v kinematic viscosity 

p density 

a surface tension 
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Subscripts : 


j 

1 

s 

V 

o 


evaluated at 
evaluated at 
evaluated at 
evaluated at 
evaluated at 


jet inlet 
liquid condition 
liquid-vapor Interface 
vapor condition 
zero Richardson number 


INTRODUCTION 

The paper is concerned with the formulation and numerical computation of 
the rate of vapor condensation on a liquid surface. The problem under consid- 
eration has direct relevance to space-based systems related to cryogenic fluid 
storage and the control of cryogenic tank pressure. The liquid-vapor interface 
location and configuration are not well defined in microgravity environments. 
Therefore, direct venting of vapor to control the tank pressure either during 
storage or during the filling process of an empty tank is not a viable option. 
Technology requirements and fluid dynamics and heat transfer problems associa- 
ted with the cryogenic fluid management in space environments are adequately 
described in numerous publications (refs. 1 to 4). In space systems, the pre- 
ferred method of storage tank pressure control is to circulate the liquid 
through a cooler, reinject the subcooled liquid into the tank and induce vapor 
condensation at the liquid surface. Therefore, the prediction of the condensa- 
tion rate in terms of the injection flow conditions and system parameters that 
may affect the condensation process is required. 

The interface condensation process which is principally controlled by the 
liquid motion beneath the Interface is of specific Interest. The liquid motion 
beneath the interface may be Induced by an axial jet of subcooled liquid. 
Recently, Sonin and his coworkers (refs. 5 to 8) have performed extensive 
experimental investigations on the condensation rate of steam on a turbulent 
liquid surface. The turbulence in the liquid side was produced by a single 
axial jet directed toward the liquid surface from below. A general empirical 
correlation for the condensation rate in terms of the liquid-side turbulence 
properties and fluid properties was developed. The effects of buoyancy on the 
condensation rate is described in terms of Richardson number. It is shown 
(refs. 6 and 8) for lower jet Reynolds numbers and a Richardson number greater 
than unity, the condensation rate is significantly reduced even for the turbu- 
lent jet. This Is due to the thermal stratification near the interface because 
of buoyancy effects. Realistically, for the ground-based condensation experi- 
ment, the buoyancy effects cannot be avoided at low jet velocities. Therefore, 
it Is essential to predict the condensation rate either analytically or numer- 
ically for the cases where the buoyancy effects are significant. 

Recently, Gerner and Tien (ref. 9) investigated the interfacial condensa- 
tion problem on a horizontal liquid surface by employing a simple axi symmetri c 
model. However, the effects of buoyancy were neglected. Lin (ref. 10) and 
Lin and Hasan (ref. 11) obtained numerical solutions of the condensation of 
vapor on a horizontal liquid surface in a cylindrical tank. The liquid motion 
beneath the interface is Induced by an axisymmetric laminar jet of subcooled 
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liquid. The effects of the jet Reynolds number, liquid Prandtl number, Jakob 
number as well as system parameters on the condensation rate are Investigated. 
Again, the effect of buoyancy was neglected. 

The present paper Is an extension of the work reported by Lin (ref. 10). 
The specific focus of the present work is to investigate the effects of buoy- 
ancy on the vapor condensation rate. 


PROBLEM FORMULATION 

The physical system and the coordinates used to analyze the problem are 
shown in figure 1. A cylindrical tank of diameter, D, is filled with liquid up 
to a height x s . The liquid surface is in contact with saturated vapor and the 
tank pressure is kept constant. Therefore, the liquid vapor interface is at a 
constant saturation temperature corresponding to the tank pressure. An axisym- 
metrlc, subcooled laminar jet is introduced at the center of the tank bottom 
and is directed toward the liquid vapor Interface. The liquid is withdrawn 
from the outer portion of the tank bottom allowing the liquid height to remain 
constant. 


In the present analysis it is assumed that the fluid flow and the temper- 
ature fields in the liquid region remain essentially unaffected by the vapor 
region (refs. 9 and 10). Therefore, the momentum equations and the energy 
equation in the liquid side can be solved simultaneously and the condensation 
rate can be determined from the resulting temperature gradient at the inter- 
face. The above simplification is based on the following assumptions: 

• The vapor is quiescent and at saturation temperature or if superheated, 
C p (T v - T s )/h f << 1. Therefore, the energy equation in the vapor phase 

need not be considered. 

• The liquid vapor interface is flat and wave free (l.e., v 2 /gD << 1 and 

u 3 /<go./p) ,/4 « l). s 

• The interface is shear free (i.e., << 1). 


The governing conservation equations in the nondimensional form can be written 
as: 


8u* 3r*v* 

3x* + r* 3r* = 


3u* 2 

3u*r*v* 

_ 3fi! + _ 

1 [s 2 u* 

a ( 

V* 9u * > i 

3x* 

+ r* 3r* 

* 3x* + B 

Re jUx* 2 + 

' 3r* V 

s 3rv 

3u*v* 

3r*v* 2 

3p* 1 

v* 1 

3 2 v* 

3 

3x* + 

r* 3r* ' 

3r* ' B Re^ r *2 T B Re^ 
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(4) 


3u*T* 3r*v*T* 1 [3 2 T* 8 

3x* + r* 3r* * B Re.. Pr[ 9x *2 + r* 8r* 



The effect of buoyancy is accounted for by retaining the body force term 
g[3(T - Tj), (Boussinesq approximation) in the x-direction momentum equation, 
where g J is the acceleration due to gravity acting normal to the surface in 
the negative x-direction. In equation (2) the effect of buoyancy is repre- 
sented by the Richardson number, Ri . 

The following boundary conditions are used to solve elliptic equations 
(1) to (4). At the centerline, the symmetric conditions are: 


v* = 0 


3u* 3T* 

3r* = 3r* = 


(5) 


At the sol id walls 


u* = v* = 0 


The bottom and side walls are assumed to be adiabatic: 


3T* 3T* 

3x* = 3r* = 


0 


( 6 ) 


(7) 


At the jet inlet 


u* = 1 v* = 0 T* = 0 


( 8 ) 


The axial velocity at the interface Is the condensation Induced velocity and is 
given by 


u 


* 

c 



B Pr Re^ 


at 



At the liquid-withdrawn plane 


u* 




(9) 


(10) 


v* = 0 



( 11 ) 


The primary parameters of physical Interest are: condensation mass flux 


m 


c 



( 12 ) 
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and the local and average condensation heat transfer coefficients and can be 
written in terms of Stanton numbers as: 


St 


c 



pUjCp B Pr Re 




8 


,0.5 



B Pr Re. 

J 


(13) 


(14) 


The nondimensional form of governing equations and the boundary conditions show 
that the average condensation Stanton number can be expressed as: 


(15) 


The effects of system parameters, Ja and Pr have been investigated in ref- 
erences 10 and 11. The focus of this study is to investigate the effects of 
buoyancy, described by the Richardson number, on the condensation process. 


( X A 

St c = f Re., Ri, Pr, Ja, jj, and 


NUMERICAL SOLUTION 

The coupled elliptic partial differential equations along with the bound- 
ary conditions are solved by a finite difference method. Calculations are 
performed with a nonuniform grid distribution with the concentration of grid 
nodes in the centerline, near wall and near Interface regions where the gradi- 
ents of the flow properties are expected to be large. The essential features 
of the numerical scheme are described In references 10 and 11. Reference 10 
shows that 72 x 41 grid nodes gives reasonably grid-independent solutions. 

The same grid distribution is used in the present computation. It should be 
noted that the coupling of momentum and energy equations is due to the buoyancy 
term as well as the condensation induced velocity, u c , at the Interface. How- 
ever, numerical results (ref. 11) show that the temperature and velocity fields 
in the liquid side, including the near-interface region, remain essentially 
unaffected by the condensation induced velocity, u c , if Ja < 0.2. Therefore, 
in the present calculation the condensation induced velocity is set equal to 
zero. The implication is that the results presented in this paper and conclu- 
sions drawn are valid only for Ja << 1 . 


RESULTS AND DISCUSSION 

Numerical calculations are performed for fixed values of the geometric 
parameters. D/d = 20 and x s /D = 1 for the jet Reynolds number ranging from 150 
to 450 and Prandtl numbers of 1.25 and 2.65. The value of the Richardson num- 
ber is varied from 0 to 0.2. 
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In the analysis, the buoyancy effect is accounted for in the Boussinesq 
approximation. The nondimensional form of the equations show that the 
Richardson number scales the buoyancy effects. In a gravitational field, buoy- 
ancy creates a warm thermal layer near the interface which reduces the conden- 
sation rate. Therefore, the velocity and temperature distributions near the 
interface and the interface temperature gradient are of Interest for the pur- 
pose of analysis. Some representative velocity and interface temperature 
gradient distributions are presented in the paper. Figures 2 and 3 show the 
radial distribution of axial velocity at two locations near the interface for 
the jet Reynolds number of 300 and Prandtl number of 1.25. The results show 
that as the Richardson number increases the jet velocity in the jet influenced 
region decreases. This occurs because the jet fluid temperature is lower than 
the surrounding fluid temperature and, therefore, in this case buoyancy opposes 
the jet flow. Figure 3 shows that for Ri = 0.2, flow reversal occurs near the 
outer part of the jet influence region. The effects of the Richardson number 
on the radial velocity at the free surface <x* = 1) and near the free surface 
(x* = 0.97) are shown In figures 4 and 5. The results show that the radial 
outflow velocity Is significantly reduced with the Increase of Richardson 
number. A reduction in the radial outflow velocity near the interface reduces 
the rate at which the heat of condensation can be transported away from the 
interface. 

For a fixed Reynolds number the radial distribution of the temperature 
gradient at the interface for different values of the Richardson number is 
shown in figure 6. The average value of the condensation heat flux is obtained 
by Integrating this profile over the entire Interface area. The results are 
shown in figures 7 and 8 for Pr = 1.25 and 2.65, respectively, where the ratio 
of the average Stanton number to Stanton number for zero Richardson number, 
St c /St co is plotted as a function of Richardson number. It can be observed 
from figures 7 and 8 that the condensation rate decreases significantly with 
Increase In Richardson number. Also, the effect of buoyancy on the condensa- 
tion rate cannot be quantified in terms of the Richardson number alone. The 
critical value of the Richardson number below which the Stanton number is not 
significantly reduced depends on Reynolds number as well as Prandtl number. 
However, analyses of the results show that for Ri <0.1, the effect of buoy- 
ancy on the vapor condensation rate can be neglected for the range of Reynolds 
numbers and Prandtl numbers covered in the present investigation. 


CONCLUSION 

The effects of buoyancy on the direct condensation of vapor on a horizon- 
tal liquid surface in a cylindrical tank are investigated numerically. The 
liquid motion beneath the interface is induced by an axisymmetric laminar jet 
of subcooled liquid. It is shown that the Richardson number scales the buoy- 
ancy effects. For fixed values of the geometric parameters, D/d = 20 and 
x s /D = 1, numerical solutions are obtained for Richardson number varying from 
0 to 0.2. 

The results show that the condensation rate is significantly reduced for 
higher values of the Richardson number. However, the reduction of the conden- 
sation rate cannot be quantified in terms of Richardson number alone since it 
also depends on both the jet Reynolds number and the Prandtl number. 
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RADIAL DISTANCE, r* 

Figure 6. - Effect of Richardson number on the temperature 
gradient profiles at the interface. 



RICHARDSON NUMBER, Ri 

Figure 7. - Average condensation Stanton number ratio as a 
function of Richardson number for Pr = 1.25. 
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Figure 8. - Average condensation Stanton number ratio as a 
function of Richardson number for Pr = 2.65. 
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